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A COMPARISON  OF  MEASURED  3.8  MICRON 
SCATTERING  FROM  NATURALLY  OCCURRING  AEROSOLS  WITH  THAT 
PREDICTED  BY  MEASURED  PARTICLE  SIZE  STATISTICS 

Introduction 

During  a recent  field  experiment  conducted  along  the  Southern  California 
coast,  four  laboratories  cooperated  in  a joint  scattering  measurements  program. 
TRW  Space  Systems,  Inc.,  Lockheed  Laboratory,  and  Avionics  Laboratory  (Air  Force) 
made  the  actual  light  measurements,  and  the  Naval  Research  Laboratory  provided 
the  particle  counting  measurements  and  the  subsequent  analysis. 

Reasons  for  the  measurement  program  were  varied,  but  the  purpose  of  this 
paper  is  to  develop  a comparison  between  actual  measured  volume  scattering 
coefficients  (at  various  angles)  and  similar  quantities  calculated  from  particle 
size  statistics,  an  assumed  index  of  refraction,  and  MIE  scattering  theory  for 
spherical  particles.  Data,  such  as  this,  taken  in  the  'middle  IR  region,  are 
rare,  and  rarer  still  is  comparison  with  predictions  derived  from  particle  size 
statistics. 

The  program  was  carried  out  during  August  and  September  1975  at  the  TRW 
Capistrano  Test  Site,  six  hundred  feet  above  sea  level.  The  weather  was 
reasonably  constant  during  the  evolution  of  a typical  day.  Each  morning  the 
sea-ii.duced  coastal  fog  would  burn  off  by  1000  hours  and  return  after  dark 
when  the  offshore  continental  breeze  subsided.  Individual  days  showed  meteoro- 
logical range  varying  from  infinity  to  a few  kilometers,  and  relative  humidities 
varying  from  30%,  (~  8 torr  H2O)  to  over  80%,  (14  torr  H2O)  with  the  average 
humidity  centered  around  50%.  (See  Table  1.)  On  several  occasions,  a hot, 
dry  wind  (the  "Santa  Anna"  condition)  blew  from  the  upper  California  desert- 
land,  elevating  temperatures  to  ~ 30°  (with  50%  or  less  relative  himidity). 

During  these  days,  numerous  brush  fires  were  also  observed  in  the  vicinity,  but 
on  only  one  occasion  were  measurements  made  with  smoke  haze  in  the  beam  path. 

The  laser  source  was  a multiline  DF  chemical  laser  with  approximately  22 
Note:  Manuscript  submitted  January  21,  1977. 
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lines  whose  wavelengths  ranged  from  3.6um  to  4.0am.  A focused  76cm  beam 
was  transmitted  over  a 1.5km  path  to  a collecting  station. 

The  various  laboratory's  instrumentation  will  be  described  elsewhere,  but 
here  it  is  sufficient  to  note  that  each  group's  field  of  view  appeared  adequate 
to  intercept  the  full  beam.  Moreover,  certain  alignment  checks  were  made  to 
insure  that  the  beam  remained  centered  in  the  FOV  of  each  device  as  the 
instrument  angle  with  respect  to  the  beam  was  changed. 

By  common  agreement  each  group  reported  its  volume  scattering  coefficient 
in  units  of  km"^  ster"^.  Overall  calibration  was  determined  from  the  instru- 
ment sensitivity  calibration  (traceable  to  a commonly  used  black-body  source), 
laser  power,  FOV,  and  path  geometry.  In  some  cases  extinction  losses  from  the 
intercepted  beam  volume  to  the  instrument  were  also  estimated. 

Since  these  issues  will  be  addressed  elsewhere,  this  paper  will  concen- 
trate on  the  intercomparison  described  earlier. 

The  Particle  Sizing  System 

The  device  used  for  aerosoL  particle  counting  was  Particle  Measuring 
Systems'  Active  Scattering  Aerosol  Spectrometer  (referred  to  as  the  "Knol lenberg" 
system  after  its  designer).  The  basic  operating  principle  is  similiar  to  other 
counters--intercepted  light  is  scattered  into  collecting  optics  and,  by  cali- 
bration with  known  size  particles,  the  height  of  the  scattered  light  pulse  is 
used  to  determine  the  particle  cross  section.  In  most  cases  the  light  source  is 
an  incoherent  lamp  brought  to  focus  in  the  sample  volume  and  blocked  by  an 
opaque  disk  just  beyond  the  point  of  focus.  A large  condenser  lens  (centered 
on  the  field  stop)  collects  scattered  light  produced  by  the  particle  passage 


2 


and  focuses  it  on  a detector.  Unless  2 u steradians  is  covered  by  the  field- 
lens-detector  geometry,  the  instrument  gives  an  apparent  particle  size  depen- 
dent on  the  particle  index  of  refraction.  In  the  usual  case,  however,  only 
the  forward  scattered  light  is  measured.  In  addition  to  a serious  refractive 
index  dependence,  signal-to-noise  considerations  limit  measurement  of  particles 
smaller  than  about  O.l^im  radius.  An  additional  practical  difficulty  with 
standard  instruments  intrudes:  the  particles  must  be  brought  into  the  focal 

volume,  necessitating  "plumbing"  with  an  attendent  modification  of  the  size 
distribution  due  both  to  wall  adherence  and  to  relative  humidity  variations. 

The  Knollenberg  device  partially  circumvents  these  difficulties  by  using 
a forced  aspiration  system  to  pass  the  particles  through  a HeNe  laser  cavity. 
The  airflow  and  cavity  optics  are  such  that  the  cavity  beam  essentially  defines 
the  sampling  volume,  thus  eliminating  small  bore  piping  and  the  attendant 
size  distribution  distortion.  Light  is  collected  and  scattered  beyond  one 
cavity  mirror  in  the  range  4°  to  22°  (170°  - 150°  backscatter)  giving  a con- 
siderably larger  collection  solid  angle  than  conventional  instruments  and 
thereby  achieving  a significant  reduction  in  refractive  index  dependence  of  the 
apparent  particle  size.  Lastly,  signal-to-noise  is  considerably  enhanced  by 
laser  flux  densities  up  to  16  times  that  for  incoherent  sources. 


Estimation  of  Scattering  Coefficient  From  Particle  Size  Distribution 

The  DF  chemical  laser  beam  was  directed  over  the  1.5  km  range  for  periods 
of  a few  seconds.  (These  "shots"  are  identified  by  a number  VLI-XXX  in  the 
subsequent  tables.)  During  an  eight  minute  interval  about  each  shot,  aerosol 
size  information  was  collected  at  a site  located  100  meter  perpendicularly 


distant  from  the  beam  and  500  meters  from  the  laser  transmitter. 


V 


A series  of  experiments  had  suggested  the  eight  minute  averaging  time. 

Long  intervals  gave  "skewed"  results  because  the  aerosol  distribution  changed 
appreciably,  while  shorter  intervals  were  unsatisfactory  because  of  the 
reduced  sample  statistics  of  the  larger  particles. 

The  range  of  the  Knollenberg  instrument  extended  from  0.10/im 
radius  to  5.0/xm  radius.  Figure  1 shows  the  arrangement  for  the  four  overlapping 
particle  size  ranges.  Pulses  from  the  counter  head  are  commutated  between  the 
four  ranges  in  a fraction  of  a second.  Pulse  height  gates,  16  per  range,  size 
the  "particles"  and  accumulate  counts  during  the  commutator  dwell.  The  results 
are  recorded  on  digital  magnetic  tape  at  one  second  intervals.  Later,  inter- 
mediate processing  averages  the  number  of  counts  in  each  "bin"  for  an  eight- 
minute  period.  At  this  stage,  numerous  numerical  experiments  were  performed 
with  the  data  in  an  effort  to  curve-fit  the  data  (on  log-log  paper)  and  to  do 
the  subsequent  numerical  integration  over  a particle  size  range  extending 
beyond  the  instrument  range.  The  idea  was  to  estimate  the  contributions  of  a 
few  large  particles  which,  because  of  their  high  forward  scattering  efficiency 
at  3.8am,  can  significantly  contribute  to  the  volume  scattering  coefficient. 

To  date,  the  approach  has  yielded  inconclusive  results.  The  uncertainties 
primarily  resulted  from  irregularities  in  the  size  spectrum  beyond  2.0am,  which 
were  inadequately  represented  by  second  or  third  order  fits.  The  best  line  of 
investigation  for  fitted  data  integration  appears  to  be  a weighted  second  order 
fit,  where  the  weights  are  proportioned  to  the  number  of  counts  in  each  size  bin. 
An  alternate  approach,  the  one  pursued  in  this  work,  eschews  the  curve- fit  and 
approximates  the  volume  scattering  coefficient  as  a "histogram  integration". 

Such  a result  is  clearly  a lower  bound  (since  the  range  of  the  counter  is  limited) 
but,  within  the  resolution  of  the  instrument,  it  reproduces  the  "bumps"  in  the 
large  particle  end  of  the  range  (which  have  a significant  effect  on  the  forward 
scattering  and  which  may  represent  quasi-mono-disperse  sea  salts  or  pollen 


particles  at  1.0  ~ 2.0iim. ) 


The  final  processing  of  the  eight  minute  average  data  consisted  of 
combining  the  particle  counts  in  each  of  the  four  ranges  to  form  720  ele- 
mentary "quantum  bins"  (each  . 005jim  wide)  extending  from  the  first  bin  thresh- 
old in  the  first  range  to  the  last  in  the  last  range.  For  each  range,  quantum 
bins  outside  of  the  range  were  filled  with  zeroes.  A composite  quantum  bin 
range  was  obtained  by  averaging  like  quantum  cells  across  the  four  ranges. 

These  quantum  cells  were  then  collapsed  by  a factor  of  12  to  form  60  equi- 
spaced  bins  approximating  the  original  instrument  resolution. 

Figures  2 and  3 show  logarithmic  plots  of  the  relative  number  density  versus 
the  relative  particle  radius  for  two  particular  experiments,  VL1-168  and  VL1-182. 
As  shown  in  the  next  section,  these  two  experiments  were  used  to  illustrate  a 
comparison  of  measured  angle  scattering  coefficients  and  predictions  based 
upon  the  particle  size  distribution  data.  The  dots  represent  "experimental" 
particle  density  measurements  processed  as  described  above,  whereas  the  solid 
lines  give  the  best  curve- fit  to  the  data  points.  There  is  about  a factor  of 
three  difference  in  the  scattering;  the  latter  shot,  VL1-182,  being  made  on  a 
clear,  hot,  relatively  dry  day.  Yet  the  shape  of  the  two  curves  is  quite 
similar,  particularly  in  the  vicinity  of  "bump"  at  about  1.8  (l.Oum  diameter) 
on  the  logarithmic  radius  scale.  Such  curves  have  been  seen  before;  the  mag- 
nitude of  the  "bump"  increases  or  decreases  (and  even  disappears  altogether 
in  some  measurements)  but  appears  to  be  real  data.  Current  speculation  suggests 
a bimodal  distribution  in  the  vicinity  of  a marine  environment  with  small  conti- 
nental haze  particles  (giving  the  initial  smooth  power-law  dependence)  super- 
imposed on  a characteristic  marine  haze  of  larger  particle  sizes  (perhaps  of 
sea  salt  plus  water  and  sea  surface  ejecta).  Evidence  of  this  type  strongly 
suggests  that  the  assumption  of  a uniform  index  of  refraction  (for  water,  r]  = 
1.363  + 0.00391  at  3.802um)  chosen  for  all  particle  sizes  is  invalid. 
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However,  no  better  data  is  available. 

Incidently,  the  appreciable  deviation  of  the  "bump"  shown  in  Figure  2 
and  3 was  one  of  the  principal  reasons  for  choosing  a histogram  integration 
technique  mentioned  above.  The  curve  fit  clearly  misses  an  appreciable  con- 
tribution to  the  scattering.  Further,  the  "graininess"  of  the  two  plots  for 
large  values  of  r is  a counter  limitation.  Even  though  the  averaging  period 
is  eight  minutes,  the  dwell  time  on  each  range  is  one  fourth  of  that.  There 
is,  therefore,  both  a count-statistics  variation  and  poor  size  resolution  due 
to  a bin  size  choice  required  to  accumulate  enough  counts. 

A table  of  scattering  efficiencies  Qsc(x>9) 

where 

x = 2jtr/X  (1) 

is  the  dimensionless  particle  size,  and  9 is  the  scattering  angle,  was 
calculated  for  the  10,000  (x,9)  pairs  in  the  range 


x [.1  (.0295)6.0]  (2) 

and 

9 [2.0°(3.5°)177°]  . (3) 


An  algorithm,  based  upon  Kerker's^^  formulation  of  the  Bessel- 

(2) 

Riccati  Series  developed  by  Hancock  and  Livingston  was  used.  Consistency 
comparison  with  RAND  tables  show  the  double  precision  calculation  is  valid  to 
4-5  significant  figures.  The  histogram  integration  is  thus  given  by 


eac<e)  = *k’3  10"3j=i  xj2  n(xj>  Q8c<xj’0)  rkia  lster~1]; 

2 3 

Xj  is  the  dimensionless  area  parameter;  n(x^)  [particles/cm  ] is  the 

number  of  counts  in  bin  j;  and  Qsc  is  the  scattering  efficiency.  The 


(4) 


wavenumber  k (=2n/X)  and  the  factor  10  scale  the  dimensions  of  the  volume 
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scattering  coefficient  to  km  ster 

It  was  found  by  measurement  that  the  laser  beam  was  unpolarized,  so  that 

the  Q is  actually  the  average  of  the  parallel  and  perpendicular  component  Q 
sc 

relative  to  the  plane  of  scattering. 

In  the  next  section  calculations  based  upon  Eq.  4 are  compared  with  the 


A Comaprlson  of  Measured  and  Estimated  Angle- 
dependent  Volume  Scattering  Coefficients 

Table  2 gives  a schematic  over-view  of  the  scattering  data  taken  from 
the  8th  through  the  19th  of  September  1975.  Thirteen  scattering  angles  were 
investigated  for  a total  of  61  measurements.  Figure  4 shows  a schematic  giving 
the  relative  placement  and  range  of  each  experimentor 1 s results.  As  shown 
later,  Avionics  Laboratory  results  appear  low  (by  a factor  of  four  or  so)  rela- 
tive to  the  other  groups  for  reasons  still  under  study.  At  this  time,  it 
does  not  appear  to  be  the  result  of  beam-receiver  field-of-view  misalignment. 
(That  was  checked  during  the  course  of  the  experiment.)  Further,  the  ratio 
to  other  groups'  data  appears  constant. 

For  the  particle  distribution  measured  during  each  of  the  shots,  the 
angular  scattering  function  was  calculated  according  to  Equation  4.  Each  re- 
sulting angular  scattering  function  estimate  was  then  summed  over  the  full 
solid  angle  to  yield  the  total  scattering  coefficient.  Then  a normalized 
angular  scattering  function  was  calculated  by  dividing  the  angular  scattering 
by  the  total  scattering  for  each  shot.  All  of  these  normalized  angular 
scattering  functions  were  then  averaged  together  at  each  angle  to  give  a 
composite  normalized  angular  scattering,  which  is  graphed  in  Figure  5.  In 
that  figure,  the  chain  dotted  and  the  chain  dashed  lines  show  the  area  bounded 
by  one  and  two  standard  deviations,  respectively,  from  the  average  value  of 
the  function  at  each  of  the  50  angles  used  in  the  cauculations.  Additionally, 
each  run's  calculated  total  scattering  was  used  to  normalize  the  scattering 
measurements  made  by  the  three  groups.  Their  normalized  measurements  are  also 
shown  in  Figure  5. 

Beyond  a direct  comparison,  two  semi- independent  quantities  (derived  from 
the  calculated  scattering  coefficient)  can  be  inferred  from  the  multiple  angle 
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data.  These  are  the  relative  shape  and  magnitude  of  the  scattering  coefficient 
curve.  The  relative  shape  can  be  inferred  by  taking  forward- to-backscatter 
ratios  and  comparing  them  to  similar  ratios  from  forward  to  backscatter  measure- 
ments. Tables  3 and  4 show  ratio  comparisons  for  5°  to  177°  and  19°  to  177°  for- 
ward-to-backscatter  angles,  respectively.  The  choice  of  angles  (which  was 
somewhat  arbitrary)  depended  primarily  on  the  quantity  of  data  available. 

Since  it  became  clear  during  the  course  of  the  experiment  that  possibly  two 
different  populations  of  particles  were  involved  (forward  scattering  from  water 
droplets  and  backscatter  from  condensation  nuclei)  the  ratios  may  prove  to  be 
a sensitive  test  of  the  uniform  refractive  index  approximation. 

A comparison  of  scattering  based  upon  the  particle  size  data  with  measure- 
ment is  essentially  statistical  in  nature  since  the  particle  sampling  instrument 
was  located  some  distance  away  from  the  sampled  region(s)  of  the  IR  beam. 

Therefore,  a significant  indicator  in  Tables  3 and  4 is  the  "ratio  of  ratios" 
average.  This  is  the  experimental  forward/backscatter  ratio  divided  by  the 
theoretical  estimate  (shownas  Re/Rt)  averaged  over  the  set  of  experiments. 

Were  agreement  perfect,  on  the  average,  the  ratio  is  unity.  In  Table  3,  the 
ratio  is  1.03.  To  interpret  this  "sample  mean"  as  drawn  an  infinite  population 
mean,  we  apply  sample  statistics.  The  result  is  that  the  hypothesis  "population  mean 


is  unity"  is  valid  with  only  a 57»  risk. 

The  same  situation  is  not  true  in  Table  4,  where  the  ratio  of  ratios  average 
is  only  0.16  and  the  hypothesis  "population  mean  is  unity"  fails. 

Therefore,  on  the  average,  the  experimental  scattering  curve  might  appear 
shallower  than  the  prediction  based  upon  particle  distribution.  Such  a result 
could  represent  the  upper  and  lower  range  limitation  of  the  counting  instrument 
(5.0Jlm)  as  well  as  an  effect  produced  by  a non-uniform  particle  index  of  refraction. 
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An  extension  of  the  "ratio  of  ratios"  to  all  other  combinations  of  angles, 
is  shown  in  Table  5.  Since  the  data  based  is  very  limited  for  some  angles,  a 
rigorous  application  of  small  sample  statistics  is  probably  not  warranted. 

Instead,  the  chief  utility  of  the  table  is  to  establish  the  internal 
consistency  of  the  various  experimental  groups.  As  alluded  to  earlier,  there 
appears  to  be  a problem  with  the  Avionics  data  as  shown  by  the  largest  deviation 
from  unity  of  the  average  shown  in  the  "Avionics-Lockheed"  and  "Avionics-TRW" 
forward/backscatter  colums.  However,  only  the  TRW-TRW  forward/backscatter 
average  admits  the  hypothesis  that  the  population  mean  is  unity  with  a 10%  risk. 

Lastly  there  appears  to  be  a trend  in  both  Tables  5 and  6 in  which  the 
ratios  tend  to  approach  unity  with  increasing  time  (experience).  If  true, 
then  some  of  the  earlier  data  showing  large  variations  could  be  partly  inter- 
preted as  systematic  error  which  was  reduced  with  increased  operator  skill. 

As  shown  previously,  there  appears  to  be  a systematic  deviation  in  the 
measured  and  calculated  angular  dependence  of  the  scattering  coefficient.  The 
next  major  question  is  how  well  the  magnitude  of  the  predicted  scattering 
coefficient  agrees  with  measurement.  The  results  for  177°  are  shown  in  Table 
6.  The  ratio  shown  suggests  that,  on  the  average,  calculation  based  upon 
particle  distribution  and  an  assumed  index  of  refraction  for  water  at  3. Sum 
is  a factor  of  two  larger  than  experiment.  The  trend  from  smaller  to  larger 
numbers  is  evident  as  time  progresses,  but  the  scatter  in  the  table  also 
increases.  Therefore,  it  is  difficult  to  judge  whether  the  trend  to  unity 
average  has  become  stable.  Runs  VL1-170,  VL1-176,  VL1-179  gave  unusually 
large  deviations  from  the  average  trend  and  were  eliminated  from  the  average. 

The  reasons  for  this  behavior  are  unclear,  but  probably  can  be  traced  to 
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in  the  input  particle  distribution  data.  Note  that  approximately  70%  of 
the  measurements  and  calculations  correlate  in  spite  of  the  different  sampling 
volumes  involved.  This  encouraging  feature  suggests  that  the  small  particle 
population  is  reasonably  uniform  in  space.  The  factor-of-two  average  disagree- 
ment between  theory  and  experiment  can  probably  be  reduced  by  co- locating 
sample  volumes,  by  obtaining  index  of  refraction  measures  for  the  small 
particles,  and  by  extending  the  range  of  the  instrument  to  larger  particle 
sizes. 

Forward  scattering  prediction  comparisons  are  shown  in  Table  7.  Here 
again  data  scatter  is  large,  and  no  trend  with  increasing  experiment  number  is 
evident.  Less  encouraging  than  the  previous  result  is  the  low  coefficient  of 
determination.  Even  though  the  predictions  are  only  33%,  low  on  the  average, 
only  28%  of  the  data  and  corresponding  calculations  correlate.  There  are  three 
possible  (and  likely)  explanations.  First,  the  Knollenberg  instrument  shows 


size  resolution  problems  beyond  2.5ixm  radius  in  which  3 . Oizm  particles  are  not 


distinguished  from  those  smaller  in  size.  Second,  the  small  number  of 


particles  in  these  size  ranges  (containing  the  more  efficient  forward  scatterers) 


leads  to  count  statistics  fluctuations.  Last,  it  may  be  that  the  large  particle 


population  is  very  inhomogeneous ly  distributed  in  space.  The  latter  issue  could 


be  resolved  by  looking  for  increased  correlation  in  a co-located  sampling 


volume  experiment 


Figure  6 is  a graph  which  shows  the  relative  contribution  whic.:  scattering 


at  a given  angle  makes  to  the  total  scattering.  The  curve  was  obtained  by 


multiplying  the  average  normalized  angular  scattering  function  by  2-'  sinO 


where  9 is  the  scattering  angle.  It  represents  the  total  scattering  at  a 


given  angle,  as  opposed  to  the  scattering  per  unit  solid  angle  which  is  graphed 
in  Fig.  5. 

Note  that  two  thirds  of  the  area  under  the  curve  lies  within  the  range 
between  10°  and  50°.  Since  the  greatest  contribution  to  the  total  scattering 
comes  from  scattering  at  these  angles,  it  is  to  be  expected  that  the  scattering 
from  this  region  will  have  greater  weight  in  determining  the  normalized  scatter- 
ing. Hence,  both  the  "necking  down"  of  the  standard  deviation  curves  in  Figure 
5 and  the  peak  in  Figure  6 are  effects  to  be  associated  with  the  solid  geometry 
of  this  problem. 

Finally,  as  can  be  seen  from  the  location  of  the  experimentally  measured 
data  points  in  Figure  5,  more  data  is  needed  in  the  range  between  10°  and  50°. 
Consideration  of  the  solid  geometry  involved  (particularly  as  evidenced  by  the 
curve  in  Figure  6)  shows  that  data  in  this  region  will  provide  a more  reliable 
measure  of  the  total  scattering. 

A Comparison  of  Experimental  Backscatter-Coefficients  with  Visible  (,65um) 

Volume  Scattering  Coefficients 

One  question  that  has  been  considered  for  some  time  is  "What  is  the 

relationship  between  visible  range  and  aerosol  scattering  in  the  IR?"  Curcio 
(4) 

and  Knestrick  considered  the  problem  in  several  papers  some  years  ago. 

Clearly,  if  the  particles  most  active  in  scattering  visible  light  also  produce 
a major  contribution  to  3.8um  backscatter  radiation  relative  to  the  (inneficient) 
backscatter  contributions  from  the  large  particles  (r  ~ 1.7um),  then  the  visible 
range  of  the  volume  scattering  coefficient  derived  from  it  ought  to  correlate  with 
3.8um  backscatter.  On  the  other  hand,  the  scattering  efficiency  of  "small 
particles"  at  3.8um  is  very  low,  so  the  question  turns  on  one  of  balance 
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"Do  many  small  particles  with  poor  scattering  efficiency  dominate  the 


backscatter  contribution  from  the  many  fewer  larger  particles?" 

Table  8 gives  the  results  and  Figure  7 shows  a plot  of  the  visible 
volume  scattering  coefficient (measured  by  optical  pyrometry)  versus  the 
Lockheed  177°  backscatter  data.  A coefficient  of  dtermination  of  65%  shows 
considerable  correlation,  suggesting  that  over  half  of  the  backscatter 
contribution  arises  from  visible-active  small  particles. 

Conclusion 

The  above  study  suggests  an  overall  picture  of  mid-IR  light  scattering 
predicted  from  haze  distributions.  The  shape  of  the  scattering  curve  re- 
produces the  experimental  measurement  ratios  to  within  a factor  of  two  or  so, 
with  a factor  of  two  overestimate  for  backscatter  coefficient  prediction  and 
with  a 38%  underestimate  for  forward  scatter  coefficients. 

It  appears  that  instrument  range  limitations  affect  both  forward  and 
backscatter  values  of  the  scattering  coefficient  with  the  index-of-refraction 
question  predominantly  open  for  the  small  particles.  Detailed  analysis  of 
the  refractive  index  variation  with  particle  size  is  required  for  further 
improvements  in  predictive  capability. 

The  amount  of  backscatter  "glory"  shown  in  Figure  5 may  be  used  to 
identify  how  irregular  the  scattering  particles  are  according  to  a recent 
paper  in  Science"*.  The  authors  claim  that  better  fits  to  angular  scatter 
data  from  microscopically  identified  "rough"  particles  follow  by  reducing 
the  amount  of  surface  wave  present  in  the  MIE  calculation.  Measurements 
given  in  this  paper  are  not  precise  enough  to  infer  particle  population  size 
regularity,  but  open  a fascinating  prospect  for  roughness  estimation  based 
upon  numerical  fitting  of  surface  wave  magnitude  as  a function  of  Bessel  index. 
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Table  3 

A Comparison  of  5°  to  177°  Forward 
to  Backscatter  Ratios 


SHOT  NO  VL  1 

re 

Rj 

re/rt 

COMMENTS 

-174 

195.4 

168.6 

1.16 

-175 

108.3 

101.6 

1.06 

-176 

527.3 

33.5 

15.7 

Theory  calculation  suspect:  eliminated 
from  average. 

-177 

-178 

182.3 

41.6 

163.8 

107.2 

1.1 

.387 

-179 

84.7 

204.7 

.42 

-18  2 

172.9 

84.4 

2.05 

1.03  AV  : .60(STD  DEVI 


rt  = fT(50)/fTn77',i 


• 1 


Table  4 

A Comparison  of  19°  to  177°  Forward 
to  Backscatter  Ratios 


SHOT  No.  VL1 

Re 

R<p 

Re/Rt 

COMMENTS 

- 

-163 

17.3 

j 

66. 9 

.259 

-164 

10.3 

80.7 

.128 

-1  ?6 

481.8 

26.8 

.056 

Theory  calculation  suspect: 

-182 

24.1 

67.7 

.356 

eliminate  from  A’vg. 

.161  A V ± .15  f STL  DEV) 

h.  - pr(19°)/8p  177° • Conclusion:  Average  forward-to-backscatter  ratio 

less  than  calculated  from  particle 
hT  * Bt  19P)/B^17Tc  distribution. 
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Table  5 

A Comparison  of  the  Various  Laboratory’s  Results 

A - RE(Angle(l)/Angle(2))/RT(Angle(l)/Angle(2)) 


Shot  No.  1 

Angle (1) /Angle  (2) 

Avionics/Lockheed 

Avlonics/TRW 

TRW/Lockheed 

TRW/ TRW 

VL1-158 

12.8/177° 

.257 

-159 

28.9/177° 

.109 

-l6o 

16°/177° 

.243 

-163 

19°/177° 

.260 

l6l°/l77° 

.204 

19/161° 

• 993 

-164 

19/177° 

.127 

161/177° 

.203 

19/161° 

.625 

-166 

101.7/177° 

.185 

-167 

101.7/177° 

.550 

-168 

l6°/177° 

.134 

101.7/163° 

• 951 

101.7/177° 

.551 

16/163° 

.231 

1630/177° 

1.05 

-170 

1677177° 

.699 

-171 

17571770 

2.13 

.-174 

57177° 

1.16 

-175 

1.07 

-177 

5°/ 177° 

1.11 

47.4/177° 

.585 

5/47.4° 

1.9 

-178 

50/1770 

.383 

-179 

5°/177° 

.420 

* 
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Avi  aics/Lockheed,  Avionics/TRW  average  ratios  show 
grt  -.est  deviation  from  unity  suggesting  internal 
inconsistency. 
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Table  6 

Comparison  of  Measured  and  Predicted  Backscatter  at  177° 


SHOT  NO. 


0 (8  )LknT1ster"'0  x 10** 


177°  Exp.  177°  Theo.  Ratio:  Exp/Theory  Comments 


Calculation  questionable 
eliminate  from  average. 


Calculation  questionable 
eliminate  from  average . 


Calculation  questionable 
eliminate  from  average. 


Coefficient  of  Determination:  70f 


.513  Average  slope 
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Table  7 

Comparison  of  Measured  and  Calculated  Forward  Scatter  at  5° 


3Ltan_1ster"1J  x 102 


SHOT  NO.  1 

5°  Exp. 

.5°  Theo. 

Ratio : Exp/Theory 

VLI-172 

3.8 

2.7 

1.41 

-173 

3-0 

2.34 

1.28 

-174 

3.L 

5.8 

.58 

-175 

1.3 

.71 

1.83 

-176 

5-t 

.0089 

651.6  Calculation  suspect:  elimin 
from  average. 

-177 

5.3 

15.5 

.374 

-176 

5.2 

1.43 

3-51 

-179 

2.5 

2.2 

• 1.13 

-lto 

£•3 

4.1 

.56 

-161 

4.9 

3.1 

1.56 

-132 

6.2 

6.2 

.756 

-183 

t .7 

5.2 

1.67 

1.33  Average  ratio 


Coefficient  of  Determination:  287! 
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Table  8 

A Comparison  of  “Visible”  Volume  Scattering  Coefficients 
(Divided  by  47r)  and  177°  Backscatter  Coefficients 


DATE 

TIME 

km' 

'3.8m  0 77°) 

km' 

' .65^/4  it 

'.65m 
RATIO  '3.8m 

1/10 

1309 

2.11-4 

8.75-3 

41.46 

1/10 

1552 

2.09-4 

9.55-3 

45.69 

1/10 

1902 

1.95-4 

9.55-3 

48.97 

0/11 

1318 

2.48-4 

1.23-2 

49.59 

1/12 

1420 

4.72-4 

2.39-2 

50.64 

0/15 

1259 

5.72-4 

4.38-2 

76.57 

1/16 

1233 

4.16-4 

8.75-3 

20.99 

1/16 

1820 

1.83-4 

7.16-3 

39.12 

1/17 

1447 

1.74-4 

5.88-3 

33.87 

1/17 

1650 

1.2-4 

5.65-3 

47.08 

1/17 

1916 

1.1-4 

3.18-3 

28.91 

1/11 

1048 

3.18-4 

7.08-3 

22.26 

1/18 

1247 

1.25-4 

4.13-3 

33.04 

1/18 

1404 

2.95-4 

3.82-3 

12.95 

0/10 

1523 

2.95-4 

7.88-3 

26.71 

■ = TRW  PREDICTIONS  FOR  MARINE  HAZE. 
WITH  21  KM  VISUAL  RANGE 
N (3.8u)=  1.37  - Olj 
N (0.55u)  = 1.53  - .006) 

= OBSERVED  RANGE  OF  VALIDATED 
SCATTERING  MEASUREMENTS 


0 20  40  60  80  100  120  140  160  180  200 

0- SCATTERING  ANGLE  (DEGREES) 

Fig.  4 — Data  range  for  each  experimental 
group  at  various  angles 
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